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Abstract. Packet-switching WDM ring networks with a hotspot transporting unicast, multicast, 
and broadcast traffic are important components of high-speed metropolitan area networks. For 
an arbitrary multicast fanout traffic model with uniform, hotspot destination, and hotspot source 
packet traffic, we analyze the maximum achievable long-run average packet throughput, which we 
refer to as multicast capacity, of bi-directional shortest-path routed WDM rings. We identify three 
segments that can experience the maximum utilization, and thus, limit the multicast capacity. We 
characterize the segment utilization probabilities through bounds and approximations, which we ver- 
ify through simulations. We discover that shortest-path routing can lead to utilization probabilities 
above one half for moderate to large portions of hotspot source multi- and broadcast traffic, and 
consequently multicast capacities of less than two simultaneous packet transmissions. We outline a 
one-copy routing strategy that guarantees a multicast capacity of at least two simultaneous packet 
transmissions for arbitrary hotspot source traffic. 

Keywords: Hotspot traffic, multicast, packet throughput, shortest path routing, spatial reuse, 
wavelength division multiplexing (WDM). 



1. Introduction 

Optical packet-switched ring wavelength division multiplexing (WDM) networks have emerged as 
a promising solution to alleviate the capacity shortage in the metropolitan area, which is commonly 
referred to as metro gap. Packet-switched ring networks, such as the Resilient Packet Ring (RPR) [H 
[2l[3], overcome many of the shortcomings of circuit-switched ring networks, such as low provisioning 
flexibility for packet data traffic [4]. In addition, the use of multiple wavelength channels in WDM 
ring networks, see e.g., [5l[6l[7l[8l[9l[10l[ITl[T2l[13], overcomes a key limitation of RPR, which was 
originally designed for a single-wavelength channel in each ring direction. In optical packet-switched 
ring networks, the destination nodes typically remove (strip) the packets destined to them from the 
ring. This destination stripping allows the destination node as well as other nodes downstream to 
utilize the wavelength channel for their own transmissions. With this so-called spatial wavelength 
reuse, multiple simultaneous transmissions can take place on any given wavelength channel. Spatial 
wavelength reuse is maximized through shortest path routing, whereby the source node sends a packet 
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in the ring direction that reaches the destination with the smallest hop distance, i.e., traversing the 
smallest number of intermediate network nodes. 

Multicast traffic is widely expected to account for a large portion of the metro area traffic due 
to multi-party communication applications, such as tele-conferences |14j . virtual private network 
interconnections, interactive distance learning, distributed games, and content distribution. These 
multi-party applications are expected to demand substantial bandwidths due to the trend to deliver 
the video component of multimedia content in the High-Definition Television (HDTV) format or in 
video formats with even higher resolutions, e.g., for digital cinema and tele-immersion applications. 
While there is at present scant quantitative information about the multicast traffic volume, there 
is ample anecdotal evidence of the emerging significance of this traffic type [T5l[16]. As a result, 
multicasting has been identified as an important service in optical networks [EJIIH] and has begun 
to attract significant attention in optical networking research as outlined in Section [1.11 

Metropolitan area networks consist typically of edge rings that interconnect several access networks 
(e.g., Ethernet Passive Optical Networks) and connect to a metro core ring [4|. The metro core ring 
interconnects several metro edge rings and connects to the wide area network. The node connecting 
a metro edge ring to the metro core ring is typically a traffic hotspot as it collects/distributes traffic 
destined to/originating from other metro edge rings or the wide area network. Similarly, the node 
connecting the metro core ring to the wide area network is typically a traffic hotspot. Examining 
the capacity of optical packet-switched ring networks for hotspot traffic is therefore very important. 

In this paper we examine the multicast capacity (maximum achievable long run average multi- 
cast packet throughput) of bidirectional WDM optical ring networks with a single hotspot for a 
general fanout traffic model comprising unicast, multicast, and broadcast traffic. We consider an 
arbitrary traffic mix composed of uniform traffic, hotspot destination traffic (from regular nodes to 
the hotspot), and hotspot source traffic (from the hotspot to regular nodes). We study the widely 
considered node architecture that allows nodes to transmit on all wavelength channels, but to receive 
only on one channel. We initially examine shortest path routing by deriving bounds and approxima- 
tions for the ring segment utilization probabilities due to uniform, hotspot destination, and hotspot 
source packet traffic. We prove that there are three ring segments (in a given ring direction) that 
govern the maximum segment utilization probability. For the clockwise direction in a network with 
nodes 1,2,...,A^ and wavelengths 1,2, ...,A (with N/A > 1), whereby node 1 receives on wave- 
length 1, node 2 on wavelength 2, . . ., node A on wavelength A, node A -|- 1 on wavelength 1, and 
so on, and with node N denoting the index of the hotspot node, the three critical segments are 
identified as: 
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(i) the segment connecting the hotspot, node to node 1 on wavelength 1, 
(a) the segment connecting node A — 1 to node A on wavelength A, and 
{in) the segment connecting node — 1 to node N on wavelength A. 

The utilization on these three segments limits the maximum achievable multicast packet throughput. 
We observe from the derived utilization probability expressions that the utilizations of the first two 
identified segments exceed 1/2 (and approach 1) for large fractions of hotspot source multi- and 
broadcast traffic, whereas the utilization of the third identified segment is always less than or equal 
to 1/2. Thus, shortest path routing achieves a long run average multicast throughput of less than 
two simultaneous packet transmissions (and approaching one simultaneous packet transmission) for 
large portions of hotspot source multi- and broadcast traffic. 

We specify one-copy routing which sends only one packet copy for hotspot source traffic, while 
uniform and hotspot destination packet traffic is still served using shortest path routing. One-copy 
routing ensures a capacity of at least two simultaneous packet transmissions for arbitrary hotspot 
source traffic, and at least approximately two simultaneous packet transmissions for arbitrary over- 
all traffic. We verify the accuracy of our bounds and approximations for the segment utilization 
probabilities, which are exact in the limit N/A oo, through comparisons with utilization proba- 
bilities obtained from discrete event simulations. We also quantify the gains in maximum achievable 
multicast throughput achieved by the one-copy routing strategy over shortest path routing through 
simulations. 

This paper is structured as follows. In the following subsection, we review related work. In 
Section [2l we introduce the detailed network and traffic models and formally define the multicast 
capacity. In Section [U we establish fundamental properties of the ring segment utilization in WDM 
packet rings with shortest path routing. In Section HI we derive bounds and approximations for the 
ring segment utilization due to uniform, hotspot destination, and hotspot source packet traffic on 
the wavelengths that the hotspot is not receiving on, i.e., wavelengths 1,2, . . . , A — 1 in the model 
outlined above. In Section[5l we derive similar utilization probability bounds and approximations for 
wavelength A that the hotspot receives on. In Section [H we prove that the three specific segments 
identified above govern the maximum segment utilization and multicast capacity in the network, and 
discuss implications for packet routing. In Section [7l we present numerical results obtained with the 
derived utilization bounds and approximations and compare with verifying simulations. We conclude 
in Section m 

1.1. Related Work. There has been increasing research interest in recent years for the wide range 
of aspects of multicast in general mesh circuit-switched WDM networks, including lightpath design. 
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see for instance [H], traffic grooming, see e.g., [20j, routing and wavelength assignment, see e.g., [2H 
[22l[23], and connection carrying capacity [24j. Similarly, multicasting in packet-switched single-hop 
star WDM networks has been intensely investigated, see for instance [25l[26l[27l[2E] ■ In contrast to 
these studies, we focus on packet-switched WDM ring networks in this paper. 

Multicasting in circuit-switched WDM rings, which are fundamentally different from the packet- 
switched networks considered in this paper, has been extensively examined in the literature. The 
scheduling of connections and cost-effective design of bidirectional WDM rings was addressed, for 
instance in [29]. Cost-effective traffic grooming approaches in WDM rings have been studied for 
instance in [SOlET]. The routing and wavelength assignment in reconfigurable bidirectional WDM 
rings with wavelength converters was examined in [32| . The wavelength assignment for multicasting 
in circuit-switched WDM ring networks has been studied in [33l[Ml[35ll36ll37l[38] . For unicast traffic, 
the throughputs achieved by different circuit-switched and packet-switched optical ring network 
architectures are compared in [39] . 




Optical packet-switched WDM ring networks have been experimentally demonstrated, see for in- 
stance [TSlBO] . and studied for unicast traffic, see for instance [5BT]l6]r7ll8]l9]ll0lllllll2lll3] . Multicasting 
in packet-switched WDM ring networks has received increasing interest in recent years [42l[T0] . The 
photonics level issues involved in multicasting over ring WDM networks are explored in pS], while a 
node architecture suitable for multicasting is studied in [44]. The general network architecture and 
MAC protocol issues arising from multicasting in packet-switched WDM ring networks are addressed 
in [40j[45]. The fairness issues arising when transmitting a mix of unicast and multicast traffic in 
a ring WDM network are examined in [l6]. The multicast capacity of packet-switched WDM ring 
networks has been examined for uniform packet traffic in [47ll48l|49l[50]. In contrast, we consider 
non-uniform traffic with a hotspot node, as it commonly arises in metro edge rings [51] . 

Studies of non-uniform traffic in optical networks have generally focused on issues arising in circuit- 
switched optical networks, see for instance [52115311541155115611571158] . A comparison of circuit-switching 
to optical burst switching network technologies, including a brief comparison for non-uniform traffic, 
was conducted in [59]. The throughput characteristics of a mesh network interconnecting routers on 
an optical ring through fiber shortcuts for non- uniform unicast traffic were examined in [60] , The 
study [M] considered the throughput characteristics of a ring network with uniform unicast traffic, 
where the nodes may adjust their send probabilities in a non-uniform manner. The multicast capacity 
of a single- wavelength packet-switched ring with non- uniform traffic was examined in [62]. In contrast 
to these works, we consider non-uniform traffic with an arbitrary fanout, which accommodates a wide 
range of unicast, multicast, and broadcast traffic mixes, in a WDM ring network. 
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Figure 2.1. Illustration of the clockwise wavelength channels of a WDM ring net- 
work with N = 1Q nodes and A = 4 wavelength channels. 

2. System Model and Notations 

We let N denote the number of network nodes, which we index sequentially by i, i = 1, . . . , A'^, in 
the clockwise direction and let := {1, . . . , N} denote the set of network nodes. For convenience, 
we label the nodes modulo N, e.g., node N is also denoted by or —N. We consider the family 
of node structures where each node can transmit on any wavelength using either one or multiple 
tunable transmitters (TTs) or an array of A fixed-tuned transmitters (FT^), and receive on one 
wavelength using a single fixed-tuned receiver {FR). 

For N = A, each node has its own home channel for reception. For N > A, each wavelength is 
shared by ry := N/A nodes, which we assume to be an integer. For 1 < ? < iV, we let 'ui denote 
the clockwise oriented ring segment connecting node z — 1 to node i. Analogously, we let Ui denote 
the counter clockwise oriented ring segment connecting node i to node i — 1. Each ring deploys the 
same set of wavelength channels {!,..., A}, one set on the clockwise ring and another set on the 
counterclockwise ring. The nodes n = X + kA with k £ {0,1, . . . ,r] — 1} share the drop wavelength 
A. We refer to the incoming edges of these nodes, i.e., the edges li^A+feA and tT^+i-f-feA, as critical 
edges on A. 

For multicast traffic, the sending node generates a copy of the multicast packet for each wavelength 
that is drop wavelength for at least one destination node. Denote by S the node that is the sender. 
We introduce the random set of destinations (fanout set) C ({1, 2, . . . ,N} \ {S}). Moreover, we 
define the set of active nodes A as the union of the sender and all destinations, i.e., A:= J^U {S}. 
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We consider a traffic model combining a portion a of uniform traffic, a portion (3 of hotspot 
destination traffic, and a portion 7 of hotspot source traffic with a, /?, 7 > and a + /? + 7 = 1: 

Uniform traffic: A given generated packet is a uniform traffic packet with probabihty a. For 
such a packet, the sending node is chosen uniformly at random amongst all network nodes 
{1,2,..., N}. Once the sender S is chosen, the number of receivers (fanout) ^ G {1, 2, . . . , N— 
1} is chosen at random according to a discrete probability distribution {lJ'i)f^^- Once the 
fanout I is chosen, the random set of destinations (fanout set) C ({1, 2, . . . , N} \ {S}) is 
chosen uniformly at random amongst all subsets of {1,2,..., N} \ {S} having cardinality I. 
We denote by Pa the probability measure associated with uniform traffic. 

Hotspot destination traffic: A given packet is a hotspot destination traffic packet with prob- 
ability (3. For a hotspot destination traffic packet, node N is always a destination. The 
sending node is chosen uniformly at random amongst the other nodes {1,2, ...,A^ — 1}. 
Once the sender S is chosen, the fanout ^ G {1, 2, . . . , — 1} is chosen at random according 
to a discrete probability distribution (f/)^^^- O^ce the fanout I is chosen, a random fanout 
subset J^' C ({1,2, . . . ,N — 1}\ {S}) is chosen uniformly at random amongst all subsets of 
{1, 2, . . . , A" - 1} \ {S} having cardinahty (/ - 1), and the fanout set is JF = :F' U {N}. We 
denote by Qjs the probability measure associated with hotspot destination traffic. 

Hotspot source traffic: A given packet is a hotspot source traffi,c packet with probability 7. 
For such a packet, the sending node is chosen to be node N . The fanout 1 < / < {N — 1) is 
chosen at random according to a discrete prob. distribution (k^)^^^. Once the fanout I is 
chosen, a random fanout set C {1, 2, . . . , A"— 1} is chosen uniformly at random amongst all 
subsets of{l,2,...,A^ — 1} having cardinality /. We denote by the probability measures 
associated with hotspot source traffic. 

While our analysis assumes that the traffic type, the source node, the fanout, and the fanout set 
are drawn independently at random, this independence assumption is not critical for the analysis. 
Our results hold also for traffic patterns with correlations, as long as the long run average segment 
utilizations are equivalent to the utilizations with the independence assumption. For instance, our 
results hold for a correlated traffic model where a given source node transmits with a probability 
p < 1 to exactly the same set of destinations as the previous packet sent by the node, and with 
probability 1 — p to an independently randomly drawn number and set of destination nodes. 

We denote by the probability measure Pa conditioned upon \T\ = I, and define and Q''^ 
analogously. We denote the set of nodes with drop wavelength A by 

(2.1) Mx:={\ + kK\ke{Q,...,'n-l}}. 
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The set of all destinations with drop wavelength A is then 
(2.2) J'x:=J'r\Mx. 

Moreover, we use the following notation: For £ G {1, . . . , — 1} we denote the probability of i 
destinations on wavelength A by /xa/, ^\,t, and ka/ for uniform, hotspot destination, and hotspot 
source traffic, respectively. Since the fanout set is chosen uniformly at random among all subsets 
of {1,2,..., A^} \ {S} having cardinality I, these usage-probabilities can be expressed by /x;, i/;, and 
K/. Depending on whether the sender is on the drop- wavelength or not, we obtain slightly different 
expressions. As will become evident shortly, it suffices to focus on the case where the sender is 
on the considered drop wavelength A, i.e., S £ M.\, since the relevant probabilities are estimated 
through comparisons with transformations (enlarged, reduced or right/left-shifted ring introduced 
in Appendix [All that put the sender in M.\. 

Through elementary combinatorial considerations we obtain the following probability distribu- 
tions: For uniform traffic, the probability for having £ G {0, . . . , / A 77} destinations on wavelength A 
is 

(2-3) ^lx,^ := > — -tk — 



i=max(l,i?) 



(7) 

For hotspot destination traffic, we obtain for wavelengths A 7^ A and ^ € {0, ...,(/ — 1) A 77} 



^-1 (V\ (N-ri-l\ 

(2-4) i^x,r.= 2^ /iv-n 

i=max(l/) V l~l ) 

as well as for wavelength A homing the hotspot and ^€ {1,...,/Ar/} 



(2.5) z^A/ := " /jv_i^ I'l- 



N-l frj-l-^ ^N-rj-^ 
1=1 \ l-l J 



Finally, for hotspot source traffic, we obtain for A 7^ A and £ € {0, ... ,1 A r]} 

N-l fri\ fN~l--ri\ 

(2-6) KA,£ := 2^ T^vZu 1^1 

/=max{l,^) \ I I 

and for A = A and £ G {0, . . . , / A (r/ - 1)} 

(2-7) KA,i := 2^ — — Ki. 

i=max(l/) V I ) 

For a given wavelength A, we denote by the probability measure Pa conditioned upon |.7-a| = i, 
and define ^ and qf^ analogously. 

Remark 2.1. Whenever it is clear which wavelength A is considered we omit the subscript A and 
write pi, qfj, or g^. 
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We introduce the set of active nodes .4a on a given drop wavelength A as 

(2.8) Ax := :Fx U {S} . 

We order the nodes in this set in increasing order of their node indices, i.e., 

(2.9) Ax = {Xx,i,X^^2, Xx,e+i}, 1 < ^^,1 < Xx,2 < ■ ■ ■ < ^a,^+i < N, 
and consider the "gaps" 

(2.10) Xx,i + (TV - X^^e+i), {Xx,2 - • • • > iXx,e+i - Xx,e), 

between successive nodes in the set Ax. We have denoted here again hy £ = £x the random number 
of destinations with drop wavelength A. 

For shortest path routing, i.e., to maximize spatial wavelength reuse, we determine the largest 
of these gaps. Since there may be a tie among the largest gaps (in which case one of the largest 
gaps is chosen uniformly at random), we denote the selected largest gap as "CLG^"" (for "Chosen 
Largest Gap"). Suppose the CLGx is between nodes Xx,i-i and Xx,i- With shortest path routing, 
the packet is then sent from the sender S to node X^^i-i, and from the sender S to node Xx^i in the 
opposite direction. Thus, the largest gap is not traversed by the packet transmission. 

Note that by symmetry, F{ui is used} = P{i?Ar is used}, and FIun is used} = F{ui is used}. 
More generally, for reasons of symmetry, it suffices to compute the utilization probabilities for the 
clockwise oriented edges. For n € {1, . . . , N}, we abbreviate 

(2.11) nx ■= is used on wavelength A. 

It will be convenient to call node N also node 0. We let Qx, = 0, .... — 1, be a random variable 
denoting the first node bordering the chosen largest gap on wavelength A, when this gap is considered 
clockwise. 

The utilization probability for the clockwise segment n on wavelength A is given by 
v 

(2.12) P (wa) = ^{a-pi {nx) ■ f^x,t + P ■ (^a) • i^x/ + 7 • ^ ("a) • nx,e) ■ 

£=0 

Our primary performance metric is the maximum packet throughout (stability limit). More specif- 
ically, we define the (effective) multicast capacity Cm as the maximum number of packets (with a 
given traffic pattern) that can be sent simultaneously in the long run, and note that Cm is given as 
the reciprocal of the largest ring segment utilization probability, i.e.. 
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3. General Properties of Segment Utilization 

First, we prove a general recursion formula for shortest path routing. 

Proposition 3.1. Let A G {1, . . . , A} be a fixed wavelength. For all nodes n € {0, . . . , — 1}, 
(3.1) F(^{n + 1)^ = P^a) +P(5 = n)-P(6?A = n). 

Proof. There are two complementary events leading to (n + 1)^: (A) the packet traverses (on wave- 
length A) both the clockwise segment 'un+i and the preceding clockwise segment 'un, i e., the sender 
is a node S ^ n, and (B) node n is the sender {S = n) and transmits the packet in the clockwise 
direction, so that it traverses segment "un+i following node n (in the clockwise direction). Formally, 



(3.2) P|^(n + l);^J = F \^nx and {n + l)^j +F = n and (n + l)^ 

Next, note that the event that the clockwise segment m„ is traversed can be decomposed into two 
complementary events, namely (a) segments Un and Un+i are traversed, and (b) segment Un is 
traversed, but not segment Un+i, i.e., 

(3.3) P (nx^ = P and (n + 1);^^ + P (^x and not (n + 1);^^ . 
Similarly, we can decompose the event of node n being the sender as 

(3.4) ¥ {S = n) = F (^S = n and {n + 1)^ + P ^5 = n and not (n + 1);^^ . 

Hence, we can express P ^(n + 1)^^^ as 

,n + l)A^ = P (^n a) - P and not (n + 1)_^ 



(3.5) +P (5 = n) - P = ra and not {n + 1) 

Now, note that there are two complementary events that result in the CLG to start at node n, such 
that clockwise segment n + 1 is inside the CLG: (i) node n is the last destination node reached by 
the clockwise transmission, i.e., segment n is used, but segment n + 1 is not used, and (ii) node n is 
the sender and transmits only a packet copy in the counter clockwise direction. Hence, 

(3.6) F{gx = n) = P (^x and not (n + 1)^ + P ^5 = n and not {n + 1)^^ . 
Therefore, we obtain the general recursion 

(3.7) F(^{n + 1)^ = F(nx) +F {S = n) -F {Qx = n) . 

□ 
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We introduce the left (counter clockwise) shift and the right (clockwise) shift of node n to be ln\^ 
and \n]^ given by 



(3.8) [n]^ :-- 



n — X 



A 



n — X 



A 



A + A. 



A + A and \n\ ^ : 

The counter clockwise shift maps a node n not homed on A onto the nearest node in the counter 
clockwise direction that is homed on A. Similarly, the clockwise shift maps a node n not homed on 
A onto the closest node in the clockwise direction that is homed on A. 

For the traffic on wavelength A, we obtain by repeated application of Proposition 13.11 

/ rv \ . . Ma~1 Ma-1 

(3.9) P (Ma)a = Ip(^a)+ E ^iS = ^)- ^iOx = ^) 



(3.10) = F\^nx)+¥{Se{n,...,\n]x-l})-^{Qx€{n,...,\n]x-l}). 

Note that the CLG on A can only start (i) at the source node, irrespective of whether it is on A, 
or (a) at a destination node on A. Consider a given node n that is not on A, then the nodes in 
{n, n + 1, . . . , \n]x — 1} are not on A. (If node n is on A, i.e., n = \ri\x, then trivially the set 
{n, n + 1,...,H.-1} is empty and F((rnl),) =Pp.).) Hence, the CLG on A can only start 
at a node in {n, n + 1, . . . , [n] a — 1} if that node is the source node, i.e., 

(3.11) FiGx € {n,...,MA-l}) =P(^A = Se {n,...,MA-l}). 

Next, note that the event that a node in {n, n+1, . . . , [n]A — 1} is the source node can be decomposed 
into the two complementary events (i) a node in {n, n + 1, ... , [n] a — 1} is the source node and the 
CLG on A starts at that node, and (ii) a node in {n, n + 1, . . . , [n] a — 1} is the source node and 
the CLG does not start at that node. Hence, 
(3.12) 

P (5 G {n, . . . , M A - 1}) = P (^A = 5 G {n, . . . , M A - 1}) + P (5 G {n, . . . , m - 1} , ^A / 5) . 
Inserting ([XTTI) and (1X121) in dSTTO]) we obtain 

(3.13) P(^([n|j^) = ¥(^nx)+r{Se{n,...,m-l},gx^S) 
which directly leads to 

Corollary 3.2. The usage of non-critical segments is smaller than the usage of critical segments, 
more precisely for n G {0, . . . , — 1}.' 

(3.14) P^a) = P(([nlj,) -P(5G{n,...,rnl,-l},(?A/5). 

To compare the expected length of the largest gap on a wavelength in the WDM ring with the 
expected length of the largest gap in the single wavelength ring, we introduce the enlarged and 
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reduced ring in Appendix El In brief, in the enlarged ring, an extra node is added on the considered 
wavelength between the A-neighbors of the source node. This enlargement results in (a) a set of ?? + l 
nodes homed on the considered wavelength, and (6) an enlarged set of active nodes Al^ containing 
the original destination nodes plus the added extra node (which in a sense represents the source 
node on the considered wavelength) for a total of ^ + 1 active nodes. The expected length of the 
largest gap on this enlarged wavelength ring with £ + 1 active nodes among r/ + l nodes homed on the 
wavelength (A) is equivalent to A times the expected length of the largest gap on a single wavelength 
ring with I = i destination nodes and one source node among N nodes homed on the ring, and (B) 
provides an upper bound on the expected length of the largest gap on the original wavelength ring 
(before the enlargement). 

In the reduced ring, the left- and right-shifted source node are merged into one node on the 
considered wavelength, resulting (a) in a set of r/ — 1 nodes homed on the considered wavelength, and 
(6) a set of ^ — 1, oi £ + 1 active nodes. The expected length of the largest gap decreases with 
increasing number of active nodes, hence we consider the case with £ + 1 active nodes for a lower 
bound. The expected length of the largest gap on the reduced wavelength ring with £ + 1 active 
nodes among rj — l nodes homed on the wavelength (A) is equivalent to A times the expected length 
of the largest gap on a single wavelength ring with I = £ destination nodes and one source node 
among nodes homed on the ring, and (B) provides a lower bound on the expected length of the 
largest gap on the original wavelength ring (before the reduction). Prom these two constructions, 
which are formally provided in Appendix [Aj we directly obtain: 

Proposition 3.3. Given that the cardinality of J-\ is £, the expected length of the CLG on wavelength 
A is bounded by: 



where g{l,N) denotes the expected length of the CLG for a single wavelength ring with N nodes, 
when the active set is chosen uniformly at random from all subsets of {1, . . . , A^} with cardinality 



The expected length of the largest gap g{l,N) f63] is given for / = 0, . . . , A^ — 1, by g{l,N) = 
Ylk=i f^'Qi,N{k), where (?/,Ar(-) denotes the distribution of the length of the largest gap. Let Pi^N{k) = 



(3.15) 



A-ff(^,r?-l) <E^(1CLGa|) < A-5(^,r/ + l) 



(/ + 1). 




k 



k-l 



(3.16) 




m=l 



m=l 
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together with the initialization <?o,Af(^) = ^N,k and qN-i,Nik) = (^i,fc, where (5jv,fc denotes the Kro- 
necker Delta. Whereby, qo,N{k) = SN,k means a ring with only one active node has only one gap 
of length A'', hence the largest gap has length A'' with probability one. Similarly, qN-i,Nik) = Si^k 
means a ring with all nodes active (broadcast case) has N gaps with length one, hence the largest 
gap has length 1 with probability one. This initialization directly implies g{0, N) = N as well as 
g{N -l,N) = 1. Obviously, we have to set g{l, N) = for I > N. 



4. Bounds on Segment Utilization for A 7^ A 

4.1. Uniform Traffic. In the setting of uniform traffic, one has for all n € {—A + A + 1, . . . , A} and 

/c € {0, . . . , r/ — 1}, for reasons of symmetry: 

(4.1) P^l^nx) =Pa(^{n + kA) 

For n G {—A + A + 1, . . . , A}, the difference between critical and non-critical edges, corresponding 
to Corollary 13.21 can be estimated by 

< Pais e{n,...,X-l},gx^S) 

(4.2) < P,(5G{n,...,A-l}) 

With shortest path routing, on average N — {\CLG\-y) segments are traversed on A to serve a 
uniform traffic packet. Equivalently, we obtain the expected number of traversed segments by sum- 
ming the utilization probabilities of the individual segments, i.e., as Yln=i (^^) ^^n=i -^o ^ 
which, due to symmetry, equals 2 Yln=i {^>^ • Hence, 

TV 

(4.3) N-Eai\CLG\^) = 2Y,Pa{nx 

n=l 

and 

N 

(4.4) Ea{\CLG\^) = N-2Y,Pa{nx) 

n=l 

(4.5) = N-2r] Pc^lkx) 
Expressing Pa ( ) using Corollary 13.21 we obtain 



(4.6) 



nx 



Eai\CLG\^) = N -2NPa(xx) +2rj Yl Pa {S e {k, . . . , \ - l},gx 7^ S) . 

^ ^ fe=-A+A+l 
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Solving for Po ( Xx] yields 



(4.7) Pa (a a) = ^ - {\CLG\^) + j Yl Pa{Se{k,...,X-l},gxj^S) 

Hence, the inequalities (14. 2p lead to 

(4-8) I - i\CLG\,) < Pa (a,) < i - ^E^ {\CLG\,) + 

Employing the bounds for Ea {\CLG\^) from Proposition 13.31 gives 



(4.9) 



1 1 X A 1 1 X A-1 



^=0 ^ ^ ' ^=0 



4.2. Hotspot Destination Traffic. The only difference to uniform traffic is that N cannot be a 
sender, since it is already a destination, i.e.. 



(4.10) 4[nx) =p^nAl5/iVj 
Using (5 = A^) = ^, we obtain 

(4.11) 4 Pa) = J^,Pi{nx)-j^/.{nx\S = N) 



(4.12) = ^/.(nA)-^g^(nA). 



Due to the factor the second term is negligible in the context of large networks. 

4.3. Hotspot Source Traffic. Since node N is the sender (and given that there is at least one 
destination node on A), it sends a packet copy over segment u„ on wavelength A if the CLG on A 
starts at a node with index n or higher. Hence, the usage probability of a segment can be computed 

as 

(4.13) q'(nx)=qUgx>n) 



for n £ {1, . . . ,N}. We notice immediately that q!^ y^xj is monotone decreasing in n. Moreover, 
for all n G {1, . . . , (r/ — 1) A + A}, Equation (|3.14p simplifies to 

(4.14) 4(nA) = g^((MA)A 

since the sender is node A'' = and consequently P (5 G {n, . . . , [n]^ — 1} , ^a 7^ 5*) = for the 
considered n G {1, ...,(?] — 1) A + A}. Since (^^) monotone decreasing in n, the maximally 
used critical segment on wavelength A is m^a. 

With node being the sender, the CLG on A can only start at the source node = 0, or at a 
destination node homed on A. If the CLG does not start at A^ = 0, the segment ux leading to the 
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first node homed on A, namely node A, is utilized. Hence, 

(4.15) ^^(^Aa) =(Z^(eA/0). 
Observe that 

(4.16) 4 {Qx = 0)<qy {Ga-X = 0) for a < |, 

which is exploited in Section I4.4[ 
Enlarging the ring leads to 

(4.17) q^^ig, = o)<q^^{g+ = o) = j^, 

since the gaps bordering node are enlarged whereas the lengths of all other gaps are unchanged. 
A right shifting of 5 yields the following lower bound: 

(4.18) q'^iGx = 0) > q'^iSr =0\X^rx)q'^{X^Tx) 



(4.19) = - — - 1 - - 

^ ^ e+l\ r] 

Thus, 



(4-20) l-^<g^(A,) < 1-^(1 



4.4. Summary of Segment Utilization Bounds and Approximation for A 7^ A. For A / A 

we obtain from ([ZT2]) and (|4T2]) 

(4.21) P [nx) = {P'^ (^a) («^a,. + ^v^M,.) + 4 {nx) (t^a,. - ^/3^A,.)) . 



£=0 

Using Corollary 13.21 for and (|4.13l) for yields 

(4.22) maxP (nx) = F ( Xa 

neA/ V / V 



i.e., the segment number A experiences the maximum utilization on wavelength A. Moreover, in- 
equality (|4.16|) yields 



rv 



(4.23) maxmaxP nA =P 1 

A^A neA/ V / > 

i.e., the first segment on wavelength 1, experiences the maximum utilization among all segments on 
all wavelengths A 7^ A. 

From (I4.2ip in conjunction with (14. 9|) and (I4.20p we obtain 



U) > ^(« + ]^/5)-^E5(Ar? + l)(a^M + ^/3^i,.) + 

^ £ / I \ 

(4-24) + E 7:^1 V""'^' " F^^'^^'V ^' 
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and 

We obtain an approximation of the segment utilization by considering the behavior of these bounds 
for large = x- Large r] imply ^±1 ~ 1 as well as ~ 1, and §{£,1] — 1) ~ (7(^, + 1)- 

Intuitively, this last relation means that the expected length of the largest gap on a ring network 
with i destination nodes among rj — 1 nodes is approximately equal to the largest gap when there 
are £ destination nodes among rj + 1 nodes. With these considerations we can simplify the bounds 
given above and obtain the approximation (valid for large rj): 

(4.26) f(i^ ~ ^{a + P)-^Y.g{e,7j)ia^ii,e + Pi^i/) + lY.7TT''^''='-P'^''- 

^ ^ ^ i=0 £=0 

5. Bounds on Segment Utilization for A = A 

For uniform traffic this case, of course, does not differ from the case A / A. 

5.1. Hotspot Destination Traffic. Since is a destination node, by symmetry it is reached by 
a clockwise transmission with probability one half, i.e., 

1 



(5.1) Q(}[Na, 2 

For hotspot destination traffic, node N can not be the sender, i.e., Q/3{S = N) = 0. Hence, by 
Proposition 13.11 

(5.2) (?A^ =1-Qp{Ga = 0). 
Moreover, we have from Corollary 13.21 with n = 1 and A = A: 

(5.3) Q/3(aa) =Qp(ia^ +Qp{S e{l,...,A-l},gA^S). 

To estimate Qp {Ga = 0), we introduce, as before, the left- resp. right-shift of 5, given by 



(5.4) L^Ja := 



A and \S~\^ :- 



A. 



Left and right shifting of S leads to the following bounds for the probability qj^ {Qa = 0), which are 
proven in Appendix IbI 

Proposition 5.1. For hotspot destination traffic, conditioning on the cardinality of J- a to he £, the 
probability that the CLG starts at node is bounded by: 
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Inserting the bounds from Proposition 15 . II and noting that < Q13 {S {1, . . . , A — 1} ,Q\ ^ S) < 

(A - l)/(2iV) leads to 

A- 1 



and 

N-l 

/ rv \ I 

(5.7) 



5.2. Hotspot Source Traffic. Since we know that is the sender and has drop wavelength A, 
we have a symmetric setting on J^a and can directly apply the results of the single wavelength 
setting [62] , 

In particular, we obtain from Section 3.1.3 in [62] for £ € {0, . . . , r/ — 1} 
(5.8) q' (na) = 



and 

(5.9) qUAA]=q'{gA^O) 



rv 



+ 1 

5.3. Summary of Segment Utilization Bounds and Approximation for A = A. Inserting the 
bounds derived in the preceding sections in (|2.12l) . we obtain 

2(r? + l) 



(5.10) f(aa) > ^a(l-i^g(A?? + l)MA/) +^/?(l- 

^ e 



and 



'^A 1/ A-l Iv^,, , \ 1^1 A-1 ^2(£r7-l) \ 

(5.11) +7 2_]— — KA,£ =:pivu, 

1=0 

whereby is given by setting A = A in (12. 3|) . Moreover, 

(5.12) r(^NA^ > ^a(^l-^Y.g{e,v + l)f^A,^ +l^=:pNl 
and 

(5.13) p(iVA^ < ^a(^l + ^^-^J2g{£,7j-l)fiA,^ +^l3=:pNu. 
Considering again these bounds for large 77, we obtain the approximations: 

(5.14) P(Aa) - -(a + /3)-£ J^5(£,r7)^A,£-/?J];— — i/A/ + 7 5^T— TKA,£=:pi« 

^ ^ ^ e=o £=i e=o 
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as well as 

(5.15) P^iVA) ~ i(a + /3)-^^5(^,r/K/=:pA^a. 

6. Evaluation of Largest Segment Utilization and Selection of Routing Strategy 

With (|4.23l) and a detailed consideration of wavelength A = A, we prove in Appendix O the main 
theoretical result: 

Theorem 6.1. The maximum segment utilization probability is 

(6.1) max max P ('nA') = max i P ( ?i ) ,P ( Aa ) ,P ( ^Va 

ne{i,...,^}Ae{i,...,A} V ; I V / V / V 

It thus remains to compute the three probabilities on the right hand side. We have no exact 
result in the most general setting (it would be possible to give recursive formulae, but these would 
be prohibitively complex). However, we have given upper and lower bounds and approximations 
in Sections 14.41 and 15.31 which match rather well in most situations, as demonstrated in the next 
section, and have the same asymptotics when r] ^ oo while A remains fixed. 

Toward assessing the considered shortest-path routing strategy, we directly observe, that P ^-^^^ 
is always less or equal to ^. On the other hand, the first two usage probabilities will, for 7 large 
enough, become larger than ^, especially for hotspot source traffic with moderate to large fanouts. 
Hence, shortest-path routing will result in a multicast capacity of less than two for large portions of 
hotspot source multi- and broadcast traffic, which may arise in content distribution, such as for IP 
TV. 

The intuitive explanation for the high utilization of the segments 1 1 and Aa with shortest-path 
routing for multi- and broadcast hotspot source traffic is a follows. Consider the transmission of 
a given hotspot source traffic packet with destinations on wavelength A homing the hotspot. If 
the packet has a single destination uniformly distributed among the other r] — 1 nodes homed on 
wavelength A, then the CLG is adjacent and to the left (i.e., in the counter clockwise sense) of the 
hotspot with probability one half. Hence, with probability one half a packet copy is sent in the 
clockwise direction, utilizing the segment Aa. With an increasing number of uniformly distributed 
destination nodes on wavelength A, it becomes less likely that the CLG is adjacent and to the left 
of the hotspot, resulting in increased utilization of segment Aa. In the extreme case of a broadcast 
destined from the hotspot to all other r/ — 1 nodes homed on A, the CLG is adjacent and to the 
left of the hotspot with probability l/rj, i.e., segment Aa is utilized with probability 1 — 1/r/. With 
probability 1—2/r] the CLG is not adjacent to the hotspot, resulting in two packet copy transmissions, 
i.e., a packet copy is sent in each ring direction. 
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For wavelength 1, the situation is subtly different due to the rotational offset of the nodes homed 
on wavelength 1 from the hotspot. That is, node 1 has a hop distance of 1 from the hotspot (in the 
clockwise direction), whereas the highest indexed node on wavelength 1, namely node {rj — 1)A + 1 
has a hop distance of A — 1 from the hotspot (in the counter clockwise direction). As for wavelength 
A, for a given packet with a single uniformly distributed destination on wavelength 1, the CLG is 
adjacent and to the left of the hotspot with probability one half, and the packet consequently utilizes 

r\- 

segment 1 1 with probability one half. With increasing number of destinations, the probability of 
the CLG being adjacent and to the left of the hotspot decreases, and the utilization of segment 1 1 
increases, similar to the case for wavelength A. For a broadcast destined to all rj nodes on wavelength 
1, the situation is different from wavelength A, in that the CLG is never adjacent to the hotspot, 
i.e., the hotspot always sends two packet copies, one in each ring direction. 

6.1. One-Copy (OC) Routing. To overcome the high utilization of the segments li and Aa due 
to hotpot source multi- and broadcast traffic, we propose one-copy (OC) routing: With one-copy 
routing, uniform traffic and hotspot destination traffic are still served using shortest path routing. 
Hotspot source traffic is served using the following counter-based policy. We define the counter Yx 
to denote the number of nodes homed on A that would need to be traversed to reach all destinations 
on A with one packet transmission in the clockwise direction (whereby the final reached destination 
node counts as a traversed node). If 1a < ''?/2, then one packet copy is sent in the clockwise direction 
to reach all destinations. If 1a > ^/2, then one packet copy is sent in the counter clockwise direction 
to reach all destinations. Ties, i.e., Yx = rj/2, are served in either clockwise or counter clockwise 
direction with probability one half. For hotspot source traffic with arbitrary traffic fanout, this 
counter-based one-copy routing ensures a maximum utilization of one half on any ring segment. 
Note that the counter-based policy considers only the nodes homed on the considered wavelength A 
to ensure that the rotational offset between the wavelength A homing the hotspot and the considered 
wavelength A does not affect the routing decisions. 

We propose the following strategy for switching between shortest path (SP) and one-copy (OC) 
routing. Shortest path routing is employed if both (I4.26P and (|5.14l) are less than one half. If (I4.26P 
or (I5.14P exceeds one half, then one-copy routing is used. For the practical implementation of this 
switching strategy, the hotspot can periodically estimate the current traffic parameters, i.e., the traffic 
portions a, /?, and 7 as well as the corresponding fanout distributions ui, and k;, Z = 1, . . . , A^— 1, 
for instance, through a combination of traffic measurements and historic traffic patterns, similar 
to [MlESlESlEHEH] . From these traffic parameter estimates, the hotspot can then evaluate (|4.26p 
and (|5l4ll . 



MULTICAST CAPACITY OF OPTICAL WDM PACKET RING FOR HOTSPOT TRAFFIC 19 

To obtain a more refined criterion for switching between shortest path routing and one-copy 
routing we proceed as follows. We characterize the maximum segment utilization with shortest path 
routing more explicitly by inserting (|4.26|) . (|5.14|) . and (|5.15|1 in (|6.2|) to obtain: 

1 a ^ 

max max F ( n\] = —(a + d) > a(f.. v^iu o 

ne{i,...,Af} Ae{i,...,A} V 2^ ' 2t] 



' e=o e=o e=i e=o J 

whereby we noted that the definition of in (12.3(1 directly implies that //a/ is independent of A. 
Clearly, the hotspot source traffic does not influence the maximum segment utilization as long as 

(6.3) 7 < 7tM,i := ^ £ ^ 

and 

1 

, 1=1 i+T^A,e 



(6.4) 7 < 7tu,A := P 



e+i' 

Thus, if 7 < 7t/ii = niiii(7t/ii,i, 'ythi,A), then all traffic is served using shortest path routing. 

We next note that Theorem 16.11 does not hold for the one-copy routing strategy. We therefore 
bound the maximum segment utilization probability with one-copy routing by observing that (|4.9p 
together with Proposition 13.21 and (|4.2p implies that asymptotically for all A € {1, . . . , A} 

1 1 

(6.5) Pa [nx) ~ - - ^E^^(^'^W- 

Hence, Pa (^^x^ is asymptotically constant. Moreover, similar as in the single wavelength case [62], 
we have 

1 



(6.6) Pf,[nx)<Pf,[NA, 2- 
Therefore, the maximum segment utilization with one-copy routing is (approximately) bounded by 

(6.7) max max P (^a) < -(a + /? + 7) - ^ V 5(^, r?)/xi,^ 

ne{l,...,Af}AG{l,...,A} V / 2 7^0 

Comparing (|6.7p with (16.21) we observe that the maximum segment utilization with one-copy 
routing is smaller than with shortest path routing if the following threshold conditions hold: 

(6.8) 7«.a = # ^""f^' 



otherwise set 74/12,1 = 00. 
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• he: 




(6.9) 




e: 



'1=1 £+i''^A,e 2 



1 ' 



otherwise set jth2,A = oo. 

If 7 > lth2 = max(7t/i2,i, 7t/i2,A), then one-copy routing is employed. 

For 7 values between •ythi and 7j/t2, the hotspot could numerically evaluate the maximum seg- 
ment utilization probability of shortest path routing with the derived approximations. The hotspot 
could also obtain the segment utilization probabilities with one-copy routing through discrete event 
simulations to determine whether shortest path routing or one-copy routing of the hotspot traffic is 
preferable for a given set of traffic parameter estimates. 



In this section we present numerical results obtained from the derived bounds and approximations 
of the utilization probabilities as well as verifying simulations. We initially simulate individual, 
stochastically independent packets generated according to the traffic model of Section [2] and routed 
according to the shortest path routing policy. We determine estimates of the utilization probabilities 
of the three segments li, Aa, and A^a and denote these probabilities by pis, pLs, and pNs. Each 
simulation is run until the 99% confidence intervals of the utilization probability estimates are less 
than 1% of the corresponding sample means. We consider a networks with A = 4 wavelength channels 
in each ring direction. 

7.1. Evaluation of Segment Utilization Probability Bounds and Approximations for 
Shortest Path Routing. We examine the accuracy of the derived bounds and approximations 
by plotting the segment utilization probabilities as a function of the number of network nodes 
A^ = 8, 12, 16, ... , 256 and comparing with the corresponding simulation results. For the first 
set of evaluations, we consider multicast traffic with fixed fanout ni = ui = ki = l/A and 
IJ'i = I'l = 1^1 = 3/(4(A^ — 2)) for / = 2,...,A^ — 1. We examine increasing portions of hotspot 
traffic by setting a = 1, /3 = 7 = for Fig. EU a = 0.6, (3 = 0.1, and 7 = 0.3 for Fig. [721 and 
a = 0.2, (3 = 0.2, and 7 = 0.6 for Fig. 17.31 We consider these scenarios with hotspot traffic dom- 
inated by hotspot source traffic, i.e., with 7 > /?, since many multicast applications involve traffic 
distribution by a hotspot, e.g., for IP TV. 

We also consider a fixed traffic mix a = 0.2, (3 = 0.2, and 7 = 0.6 for increasing fanout. We 
consider unicast (UC) traffic with ;Ui = z/i = ki = 1 in Fig. 17.41 mixed traffic (MI) with = vi = 
Ki = 1/2 and ^ = ui = m = 1/{2{N - 2)) for / = 2, . . . , A^ - 1 in Fig. EH muhicast (MC) traffic 



7. Numerical and Simulation Results 
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Figure 7.1. Segment utilization probability as a function of number of Nodes N for 
a = 1, /3 = 0, 7 = 0, and fii = ui = ki = 1/4 and jii = vi = ni = 3/(4(A^ — 2)) for 
Z = 2, ...,iV- 1. 



0.6 

■H 

^ 0.5 

O 

U 0.4 

C 

o 0.3 

■H 

S 0.2 

■H 

0.1 



plu 
pis 
pia 
pll 



50 100 150 200 
Number of Nodes 



ll 



. 5 

. 45 

. 4 
.35 

. 3 
.25 

0.2 
.15 

. 1 
.05 





















































pLu 














pLs 








pLa 
nil 























































50 100 150 200 
Number of Nodes 



(b) P 



0.4 

0.35 
0.3 

.25 
0.2 

.15 
. 1 

. 05 

















































































pNu 

pNs 








pNa 








pNl 



50 100 150 200 
Number of Nodes 



(c) P 644 



Figure 7.2. Segment utilization probability as a function of number of Nodes for 
a = 0.6, (3 = 0.1, 7 = 0.3, and ^1 = 2^1 = ^1 = 1/4 and = vi = m = 3/(4(7V - 2)) 
for / = 2,...,iV- 1. 
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Figure 7.3. Segment utilization probability as a function of number of Nodes N for 
a = 0.2, 13 = 0.2, 7 = 0.6, and = vi = ki = 1/4 and ^ = vi = m = 3/(4(7V - 2)) 
for Z = 2, ...,iV- 1. 



with ^jLi = vi = Ki = 1 / {N — 1) iox I = 1, . . . , N — 1 m Fig. 17.61 and broadcast (BC) traffic with 
IJLN-i = I'N-i = KN-1 = 1 in Fig. [721 
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Figure 7.4. Segment utilization probability as a function of number of Nodes N for 
a = 0.2, P = 0.2, 7 = 0.6, and unicast (UC) traffic with /xi = i^i = ki = 1. 
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Figure 7.5. Segment utilization probability as a function of number of Nodes N for 
a = 0.2, (3 = 0.2, 7 = 0.6, for mixed (MI) traffic with fxi = i^i = ki = 1/2 and 
Hl = i^l = Ki = l/(2(iV - 2)) for i = 2, . . . , - 1. 
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Figure 7.6. Segment utilization probability as a function of number of Nodes A'' for 
a = 0.2, P = 0.2, 7 = 0.6, for multicast (MC) traffic with ^ = i^i = ki = 1/{N - 1) 
for / = 1,...,A^-1. 



We observe from these figures that the bounds get tight for moderate to large numbers of nodes 
N and that the approximations characterize the actual utilization probabilities fairly accurately for 
the full range of N. For instance, for N = nodes, the difference between the upper and lower 
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Figure 7.7. Segment utilization probability as a function of number of Nodes N for 
a = 0.2, P = 0.2, 7 = 0.6, for broadcast (BC) traffic with ij,n-i = ^N-i = ^^N-i = 1- 
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bound is less than 0.06, for = 128 this difference shrinks to less than 0.026. The magnitudes 
of the differences between the utilization probabilities obtained with the analytical approximations 
and the actual simulated utilization probabilities are less than 0.035 for = 64 nodes and less than 
0.019 for N = 128 for the wide range of scenarios considered in Figs. 17.11 - 17. 71 (When excluding the 
broadcast case considered in Fig. 17.71 these magnitude differences shrink to 0.02 for = 64 nodes 
and 0.01 for = 128 nodes.) 

For some scenarios we observe for small number of nodes A^ slight oscillations of the actual utiliza- 
tion probabilities obtained through simulations, e.g., in Fig. I7.4r a) and 17.51 (a). More specifically, we 
observe peaks of the utilization probabilities for odd r] and valleys for even rj. These oscillations are 
due to the discrete variations in the number of destination nodes leading to segment traversals. For 
instance, for the hotspot source unicast traffic that accounts for a 7 = 0.6 portion of the traffic in 
Fig. 17.4( a). the utilization of segment 1 1 is as follows. For even rj, there are r]/2 possible destination 
nodes that result in traversal of segment 1 1, each of these destination nodes occurs with probability 
1/{N — 1); hence, segment 1 1 is traversed with probability A^/[2A(A^ — 1)]. On the other hand, for 
odd 7], there are {rj + l)/2 possible destination nodes that result in traversal of segment li; hence, 
segment 1 1 is traversed with probability (A^ + A)/[2A(A^ — 1)]. 

Overall, we observe from Fig 17.11 that for uniform traffic, the three segments governing the max- 
imum utilization probability are evenly loaded. With increasing fractions of non-uniform traffic 
(with hotspot source traffic dominating over hotspot destination traffic), the segments li and 44 
experience increasing utilization probabilities compared to segment 644, as observed in Figs. 17.21 
and 17.31 Similarly, for the non-uniform traffic scenarios with dominating hotspot source traffic, we 
observe from Figs. I7.4ll7.7l increasing utilization probabilities for the segments li and 44 compared 
to segment 644 with increasing fanout. (In scenarios with dominating hotspot destination traffic, not 
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Figure 7.8. Maximum segment utilization probability as a function of fraction of 
hotspot source traffic 7 (with a = 1 — /3 — 7) for shortest path (SP) and one-copy 
routing (OC) for fixed fraction of hotspot traffic /3 for unicast (UC) traffic, mixed 
(MI) traffic, multicast (MC) traffic, and broadcast (BC) traffic. 

shown here due to space constraints, the utilization of segment 644 increases compared to segments 
1 1 and 44.) 

In Figs. 17.31 17.61 and 17.71 the utilization probabilities for segments li and 44 exceed one half for 
scenarios with moderate to large numbers of nodes (and correspondingly large fanouts), indicating 
the potential increase in multicast capacity by employing one-copy routing. 



7.2. Comparison of Segment Utilization Probabilities for SP and OC Routing. In Fig. 17.81 
we compare shortest path routing (SP) with one-copy routing (OC) for unicast (UC) traffic, mixed 
(MI) traffic, multicast (MC) traffic, and broadcast (BC) traffic with the fanout distributions defined 
above for a network with N = 128 nodes. The corresponding thresholds 'jthi and 'yth2 are reported 
in Table [TJ For SP routing, we plot the maximum segment utilization probability obtained from the 
analytical approximations. For OC routing, we estimate the utilization probabilities of all segments 
in the network through simulations and then search for the largest segment utilization probability. 

Focusing initially on unicast traffic, we observe that both SP and OC routing attain the same 
maximum utilization probabilities. This is to be expected since the routing behaviors of SP and 
OC are identical when there is a single destination on a wavelength. For /? = 0.1, we observe with 
increasing portion of hotspot source traffic 7 an initial decrease, a minimum value, and subsequent 
increase of the maximum utilization probability. The value of the maximum utilization probability 
for 7 = is due to the uniform and hotspot destination traffic heavily loading segment 644. With 
increasing 7 and consequently decreasing a, the load on segment 644 diminishes, while the load on 
segments li and 44 increases. For approximately 7 = 0.4, the three segments li, 44, and 644 are 
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Table 1 . Thresholds ^thi and 7^/12 for scenarios considered in Fig. 17.81 

about equally loaded. As 7 increases further, the segments li and 44 experience roughly the same, 
increasing load. For (3 = 0.2 we observe only the decrease of the maximum utilization probability, 
which is due to the load on segment 644 dominating the maximum segment utilization. For this 
larger fraction of hotspot destination traffic we do not reach the regime where segments 1 1 and 44 
govern the maximum segment utilization. 

Turning to broadcast traffic, we observe that SP routing gives higher maximum utilization proba- 
bilities than OC routing for essentially the entire range of 7, reaching utilization probabilities around 
0.9 for high proportions of hotspot source traffic. This is due to the high loading of segments li 
and 44. In contrast, with OC routing, the maximum segment utilization stays close to 0.5, resulting 
in significantly increased capacity. The slight excursions of the maximum OC segment utilization 
probability above 1/2 are due to uniform traffic. The segment utilization probability with uniform 
traffic is approximated (not bounded) by (16. 5p . making excursions above 1/2 possible even though 
hotspot destination and hotspot source traffic result in utilization probabilities less than (or equal) 
to 1/2. 

For mixed and multicast traffic, we observe for increasing 7 an initial decrease, minimum value, and 
subsequent increase of the maximum utilization probability for both SP and OC routing. Similarly 
to the case of unicast traffic, these dynamics are caused by initially dominating loading of segment 
644, then a decrease of the loading of segment 644 while the loads on segments 1 1 and 44 increase. 
We observe for the mixed and multicast traffic scenarios with the same fanout for all three traffic 
types considered in Fig 17.81 that SP routing and OC routing give essentially the same maximum 
segment utilization for small 7 up to a "knee point" in the SP curves. For larger 7, OC routing gives 
significantly smaller maximum segment utilizations. We observe from Table [T] that for relatively 
large fanouts (MC and BC), the ranges between ^thi and 'yth2 are relatively small, limiting the need 
for resorting to numerical evaluation and simulation for determining whether to employ SP or OC 
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Figure 7.9. Maximum segment utilization probability as a function of fraction of 
hotspot source traffic 7. Fixed parameters: N = 128 nodes, (3 = 0.4, /i; = 1/16 for 
1 = 1,..., 16. 
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Table 2. Thresholds ^thi and ^th2 for scenarios considered in Fig. 17.91 



routing. For small fanouts (UC and MI), the 7 thresholds are far apart; further refined decision 
criteria for routing with SP or OC are therefore an important direction for future research. 

We compare shortest path (SP) and one-copy (OC) routing for scenarios with different fanout 
distribution for the different traffic types in Fig. 17.91 for a ring with = 128 nodes. We observe 
from Fig. I7.9r a) that for hotspot source traffic with large fanout, SP routing achieves significantly 
smaller maximum segment utilizations than OC routing for 7 values up to a cross-over point, which 
lies between 74^1 and 7th2- Similarly, we observe from Fig. l7.9r b) that for small 7, SP routing achieves 
significantly smaller maximum segment utilizations than OC routing for hotspot destination traffic 
with small fanout. For example, for unicast hotspot destination traffic (i.e., ui = 1), for 7 = 0.21, 
SP routing gives a multicast capacity of Cm = 3.72 compared to Cm = 3.19 with OC routing. By 
switching from SP routing to OC routing when the fraction of hotspot source traffic 7 exceeds 0.31, 
the smaller maximum utilization probability, i.e., higher multicast capacity can be achieved across 
the range of fractions of hotspot source traffic 7. 
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8. Conclusion 

We have analytically characterized the segment utilization probabilities in a bi-directional WDM 
packet ring network with a single hotspot. We have considered arbitrary mixes of unicast, multicast, 
and broadcast traffic in combination with an arbitrary mix of uniform, hotspot destination, and 
hotspot source traffic. For shortest-path routing, we found that there are three segments that can 
attain the maximum utilization, which in turn limits the maximum achievable long-run average 
multicast packet throughput (multicast capacity). Through verifying simulations, we found that our 
bounds and approximations of the segment utilization probabilities, which are exact in the limit for 
many nodes in a network with a fixed number of wavelength channels, are fairly accurate for networks 
with on the order of ten nodes receiving on a wavelength. Importantly, we observed from our segment 
utilization analysis that shortest-path routing does not maximize the achievable multicast packet 
throughput when there is a significant portion of multi- or broadcast traffic emanating from the 
hotspot, as arises with multimedia distribution, such as IP TV networks. We proposed a one-copy 
routing strategy with an achievable long run average multicast packet throughout of about two 
simultaneous packet transmissions for such distribution scenarios. 

This study focused on the maximum achievable multicast packet throughput, but did not con- 
sider packet delay. A thorough study of the packet delay in WDM ring networks with a hotspot 
transporting multicast traffic is an important direction for future research. 



Appendix A. Definition of Enlarged and Reduced Ring as well as of Left (.4^) 

AND Right Shifting (.4^) of Set of Active Nodes 

In this appendix, we first define the enlarging and reducing of the set of "A-active nodes" Ax := 
J^x U {S}. Suppose that \J^x\ = £. Depending on the setting, and with Aix denoting the set of nodes 
homed on a given wavelength A, the set J-x is chosen uniformly at random among 

• all subsets of A^a (uniform traffic and for A 7^ A also hotspot destination and source traffic), 
or 

• all subsets of Aix that contain N (hotspot destination traffic for A = A since N is always a 
destination for hotspot destination traffic), or 

• all subsets of Mx that do not contain N (hotspot source traffic for A = A since N is always 
the source for hotspot source traffic). 

Assuming S ^ Mx, we define: 
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enlarged ring: We enlarge the set M\ by injecting an extra node homed on A between \_S\^ 
and [5];^ (and correspondingly A — 1 nodes homed on the other wavelengths). After a re- 
numeration starting with at the new node (which is accordingly homed on wavelength 
A after the re-nmneration) , we obtain M.K,r)+i '■= {mA [ m G {0, ...,r/}}. We define the 
enlarged set to equal the renumbered set T\ united with the new node. This procedure 
leads to a random set of active nodes A^^ = J^^ that is uniformly distributed among all 
subsets of MK,r)+i with cardinality {I + 1) containing node 0. Note that the largest gap of 
the enlarged set is larger or equal to the largest gap of ^a- 



enlarge 




Figure A.l. Example of enlarging for N = 16, A = 4. The sender homed 
on wavelength 1 is represented by S in the left illustration. The nodes of A^3 are 
indicated by longer tick marks and the nodes of JF3 are circled. The enlarged ring has 
a total of + A = 20 nodes, with 77 + 1 = 5 nodes homed on each wavelength. The 
added node on wavelength 3 is numbered with and lies between the former [5"]^^ 
and \S'\^. 

reduced ring: We transform the set M.\ by merging the nodes \_S\^ and [S']_)^ to a single active 
node (eliminating the A — 1 nodes inbetween). After re-numeration starting with at this 
merged node, we obtain an active set on M^^r)-!- 

Depending on the cardinality of fl {[-^J^;^ , r'S']^;^} the new active set has I, or 

i — \ elements. More specifically, if neither the left- nor the right-shifted source node was 
a destination node, then |^^| = i + 1. If either the left- or the right-shifted source node 
was a destination node, then = L If both the left- and right-shifted source node were 
destination nodes, then |^^[ = I — \. In each of these cases A^^ is uniformly distributed 
among all subsets of AlA,r?-i with cardinality \A^\ that contains node 0. 
Observe that in all cases, the largest gap of A~^ is smaller or equal to the largest gap of A\. 

We also define the following transformations: 
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reduce 




Figure A. 2. Example of reducing for iV = 16, A = 4. The sender is represented by 
S and the nodes of Ms have longer tick marks. The nodes of J^s are circled. The 
nodes IS\^ and \S~\^ (as well as the 3 nodes inbetween) are merged into the node 
numbered in the right illustration. 

Left (counter clockwise) shifting: Since S is uniformly distributed on {1, . . . ,N}, the set 

(A.l) A'^:=TxU{[S\J 

is a random subset of M\. We can think of A'^ as being chosen uniformly at random among 
all subsets of M\ having cardinality {A^l and subject to the same conditions as J^\. 
Notice that \A^\ = \ J^\\ if \_S\^ G J^x and |^^| = \ + 1 otherwise. 



left shift 




Figure A. 3. Example of left shifting for = 16, A = 4. The destination nodes 
are circled on the left, and the active nodes are circled on the right. The nodes are 
renumbered after the shifting, starting with the former sender at 0. Also, the active 
nodes is renumbered, starting with Xi > 0, the first active node after the former 
sender. The former sender is therefore the last active node, i.e., A4 = 0. 

Right (clockwise) shifting: Analogously we define 

(A.2) Ax:=TxU{\S]^}. 

This is a random set chosen uniformly at random among all subsets of Mx having cardinality 
|v43^| and subject to the same conditions as J^x- 
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right shift 




Figure A. 4. Example of right shifting for = 16, A = 4. After renumbering, the 
former sender is A3 = 0. 



Appendix B. Proof of Proposition 15.11 on Bounds for Probability that CLG starts 

AT Node for Hotspot Destination Traffic for A = A 

Proof. Conditioned on S" € A^a, we obtain 



(B.l) 



(^A = 1 S e Mk) 



Hence, we only have to consider the case S ^ A^a- We will not explicitly write down this condition. 

Consider the right shifting and denote by the starting point of the chosen largest gap of . 
Since A^ = is the only fixed active node, the first gap, i.e., {0, . . . , Aa,i}, is the only one that never 
shrinks, while the last gap, i.e., {Aa/+i, . . . , A^}, is the only one that never grows. Therefore, 

(B.2) 4(aA = 0) < q'p{g^=Q). 

For reasons of symmetry, we have 

(B.3) 
and 



+ 1' 



(B.4) 



The remaining probabilities can be computed as ([5'] a ^ -^a I S ^ AIa) = ^, leading to the desired 
upper bound. 

Analogously, the left shifting yields a lower bound, namely 



(B.5) 



'?^(^a = o| L5Ja/o) > (?^(ar = o| L^Ja/o). 
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Again for reasons of symmetry, we obtain 

(B.6) 9^(^r = 0| [S\^iT^) = j^ 

and 

(B.7) 4(^r = 0| LSJaG^a\{0}) = ^. 

Finally, we have, of course, g;^ ( [5J^ G J^a | S" ^ Mk) = | and qp{lS\f^ G J^a \ {0} \ S ^ Ma) = 
I -I. □ 



Appendix C. Proof of Theorem 16.11 on the Maximal Segment Utilization 

Proof. Due to Equation (|4.23p . we only have to prove the case of drop wavelength A. 

Corollary 13.21 tell us that it suffices to consider the critical segments. Let n = 5A with 1 < 5 < rj 
be a critical segment for A. Analogously to the proof of the domination principle in |62| . we reduce 
the domination principle for hotspot destination traffic to the statement 

(C.l) g^(n>gA>n-A)>-^4(e?A>n-A), 

and for hotspot source traffic to: 

(C.2) qi. (Ga = n)> -^a' {Qa > n) . 



Figure C.l. Illustration of statement (C.l): the mean slope of a certain period is 
bigger or equal than the mean slope over all later periods 

In the 7 (hotspot source traffic) setting, we know that N is the sender, and thus Aa C M.a- 
Hence, we do not need to consider the nodes on the other drop wavelengths and the proof is exactly 
the same as in the single wavelength case [62], see also figure C.2. 

We will now use the same strategy for the more complicated proof in the /3 (hotspot destination 
traffic) setting. Let Kn denote the number of active nodes finding themselves between the nodes 
and n (clockwise), i.e., 

(C.3) Kn:=\AAn{l,...,n-A}\. 
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Figure C.2. Gamma setting: the usage probability stays constant on non critical edges 

For k S {0, . . . , (n — 1) A — 1)} we denote q^'^ for the probability measure conditioned on 
Kn = k. We denote again n = for (5 G {1, . . . , 77 — 1}. We will show that 

(C.4) q'p{n-A<gA<n)> (^a > n - A) . 




Figure C.3. Beta setting: the usage probability changes along each segment 

In case that 5 G M.\ we can again use the proof of the one wavelength scenario. This is also true 
if 5 G {1, . . . , n — A}, since we do not claim anything about these nodes. Hence, we only have to 
investigate the case 5'G{n — A + 1,..., N} \ Ma- From now on we assume this to be the case. 

We decompose the left hand side into two parts, 

(C.5) 4 (n - A < < n) = (e?A = n) + (6^A = 5, n - A < 5 < n) . 

For the first summand of (IC.5I) . we proceed similarly to the case of a single wavelength, namely 

qp'' {Qa =n) = q^ {Qa = n,QA> n, [5]^^^ 7^ n, n G J^a) 

= qp (^A = 1 > [5Ja ^n, n e JTa) x 
xq'j^'' {gA>'n \ [S\j^ j^n,ne Fa) 
(C.6) xqf{\_S\^^n,nerA). 
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We obtain 

qf{QK = n\G>n,[S\^^n,n(^J^^) = (^^ = 1 [ [5J^ ^ 1, 1 G ^a) 

(C.7) > 4:^_„+A (^r = 1 1 l^Ja 7^ 1, 1 € ^a) . 

This probability can be computed precisely 

4:^-n+A (^r = 1 1 l5Ja^i,ig-^a) 

= q'-fn (Sa = 0) 4^-„+A (L^Ja ^ I L^Ja / 1, 1 G -^a) + 

(Sa = 0) 4:^_„+A (L5Ja e ^A I lS\^^l,le ^a) 

1 ^-^-^V 1 £-k-l 



i-k \ r]-6 J £-k-l i] - 6 



(C-8) =T^fl+ ^ 



i — k \ V ~ ^ , 

We now use the fact that, conditionally on 5 S {n — A + 1, . . . , A^} \ ^A\^ 
(C.9) q'p'ilSlA ^neTA) = 4' (L^Ja / n) gj'' (n G ^a) 



and 



(C.IO) gJ''(eA>n| L^J^y^nG^A) 



qp {QA>n\ne J^a) 



qp''ilS\A¥=n) 

ri~5 



Hence, we obtain with q^^ {n G ^a) = ^ ri~s^ ^^^^ 



4''' i^A = n) > ^^-3^9^''' (^A > n I n G J^a) 



7] — 6 

For the second part of (jC.Sp . we obtain 

q^l^^ {Qa e Is\n) = (C/a = S,Qi^> S, \S]^ = n, n ^ J^a) 

= q^f^'' {Ga = S\gA> S, \S]jy = n,n^ Ta) x 
(C.12) xq'/ {gA>S\ \S]^ = niTA) qt{\S^A = n,ni Ta) 

We have 

q^p" {Qa = S\Qa> S, \S]^ = n, n ^ Ta) 
(C.13) = , A (^A = ^ I \S^A = 1, 1 ^ -^a) . 
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Now, we use that ~ l-^A + 1| for l"/?]^ ^ Tt^. Hence, we obtain 



> 4:i^-n+A (^r = 1 1 r^iA = i, i^^a) 



Note that, conditioned on S E {n — A + 1, . . . , N} \ Ma, we have 

(C.15) Qf{\SU = niTK) = 4'' (r^lA = n) q'^' (n ^ J^a) 

and 



qf{GA>S\ \S]j, = n) 



(C.16) gJ'*^(^A>5| \S]^ = n^J^A) 



Summarizing, we obtain, using q^ (["^Ia = = that 

(C.17) gj'^ (^a ^h\n) > (^a ^ ^ I T^Ia = n) . 

It remains to show that 

XQp'' iGA>n\ne J^a) + 

1 / e-k-1 
1- 



£ — k \ rj — 6 

(C.18) xqf {gA>S \ \S^^ = n, n ^ ^a) . 



This can be shown by 



q f (^A > n - A) 
= Q^{QK>i^Xk+i=iK)qf{Xk+i=iK) + 



i=S 
A-1 



+ Yl it (^^A > iA - A I =ik-\) qf (Xfc+i =ik-\) 

A=l 

-(-*■" ^t (^A > n I Xfc+i = n) + 

(C.19) +_l_<^J.'=(g^>S| r5l^ = n). 



For the last inequahty, we used that for i G {5, . . . , r/ — 1} and A G {0, . . . , A — 1} 

qf {Qa>iA-\\ Xk+i =iA-X) 
(C.20) < gj'*^ {gA>n-X\ Xk+i = n-X) 
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and, for reasons of symmetry, 

(C.21) gJ.'=(X,+iG^A) = l-^. 

The last step we need is a comparison of (IC.ISP and (|C.19I) . The only difference arises, when both 
of the events, {n € ^a} and {[S]^ = n}, take place. Then, 

(C.22) q1^' igA>S\\S]^ = ne Ta) = q'^ {QA>n\\S^^=n(^ ^a) • 

This occurs with probability g^'^ (n G J^a \ ["5*1 a = ^) = ^'^^ explains the additional factor in 

the decomposition (IC.lSp . □ 
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